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We demonstrate a spatially resolved autocorrelation mea- 
surement with a Bose-Einstein condensate (BEC) and mea- 
sure the evolution of the spatial profile of its quantum me- 
chanical phase. Upon release of the BEC from the mag- 
netic trap, its phase develops a form that we measure to be 
quadratic in the spatial coordinate. Our experiments also re- 
veal the effects of the repulsive interaction between two over- 
lapping BEC wavepackets and we measure the small momen- 
tum they impart to each other. 



A trapped Bose-Einstein condensate has unique 
value as a source for atom lasers 0] and matter-wave in- 
terferometry [|| because its atoms occupy the same quan- 
tum state, with uniform spatial phase. However, when 
released from the trapping potential, a BEC with repul- 
sive atom-atom interactions expands, developing a non- 
uniform phase profile. Understanding this phase evolu- 
tion will be important for applications of coherent matter 
waves. We have developed a new interferometric tech- 
nique using spatially resolved autocorrelation to measure 
the functional form and time evolution of the phase of 
a BEC wavepacket expanding under the influence of its 
mean field repulsion. 

In 1997, the coherence of weakly interacting BECs 
was demonstrated by releasing two spatially separated 
condensates and observing their interference Sub- 
sequent experiments have further investigated conden- 
sate coherence properties. One || used velocity-resolved 
Bragg diffraction || to probe the momentum spectrum 
of trapped and released BECs. A complementary ex- 
periment ]7j that used matter-wave interferometry can 
be interpreted as a measurement of the spatial correla- 
tion function, whose Fourier transform is the momentum 
spectrum. These experiments showed that a trapped con- 
densate has a uniform phase, and a released condensate 
develops a non-uniform phase profile. (Recently the influ- 
ence of non-zero temperature on coherence properties was 
also investigated Q). The experiments reported in this 
Letter combine spatial resolution and interferometry to 
measure the functional form of the time-dependent phase 
profile of a released condensate. We also make the first 
measurement of the velocity imparted to two equal BEC 
wavepackets from their mutual mean- field repulsion ||] . 

We perform our experiments with a condensate of 
1.8(4) x 10 6 @ sodium atoms in the 35 1/2 , F = 1, 



nip = —1 state. The sample has no discernable non- 
condensed (i.e. thermal) component. The condensate is 
prepared following the method of Ref. || and is held in 
a magnetic trap with trapping frequencies uo x = y/2ui y = 
2uj z = 27rx27 Hz. Using a scattering length of a = 2.8 
nm, the calculated Thomas-Fermi diameters [|ll| are 47 
/Ltm, 66 /jm, and 94 /im, respectively. 

We release the BEC from the magnetic trap and it 
expands, driven mostly by the mean-field repulsion of 
the atoms. This expansion implies the development of 
a nonuniform spatial phase profile (recall that the veloc- 
ity field is proportional to the gradient of the quantum 
phase). After an expansion time To, we probe the phase 
profile with matter- wave Bragg interferometry [Ea-jLj). 
Our interferometer splits the BEC into two wavepackets 
and recombines them with a chosen overlap, producing 
interference fringes, which we measure with absorption 
imaging From the dependence of the fringe spac- 

ing on the overlap, we extract the phase profile of the 
wavepackets. 
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FIG. 1. Space-time diagram of the experiment. Three 
optically induced Bragg-diffraction pulses form the interfer- 
ometer. The condensate is released for a time To before the 
first Bragg pulse. The centers of ipA and t/>b are separated 
by Sx at the time of the third Bragg pulse, which splits them 
into ipAl, ipBi, and tpA2, ipB2- Before imaging the atoms, we 
allow the output ports to separate for a time Tz ~ 2 ms. The 
image shows the output ports when To = 3 ms, T\ = 1 ms, 
and T2 = 1.3 ms. 

Our atom interferometer Jl4|] consists of three 
optically-induced Bragg-diffraction pulses applied succes- 
sively in time (Fig. 1). Each pulse consists of two counter- 
propagating laser beams whose frequencies differ by 100 
kHz. They are detuned by about —2 GHz from atomic 
resonance (A = 2-7r/fc = 589 nm) so that spontaneous 
emission is negligible. The first pulse has a duration 
of 6 fj,s and intensity sufficient to provide a 7r/2 pulse, 
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which coherently splits the BEC into two wavepackets, 
ipA and tpB- The wavepackets have about the same num- 
ber of atoms and only differ in their momenta: p = and 
p = 2hk. At a time Ti = 1 ms after the first Bragg pulse, 
the two wavepackets are completely separated and a sec- 
ond Bragg pulse (a tt pulse) of 12 us duration transfers 
TpB to a state with p w and ipA to p ~ 2hk [jl6). Af- 
ter a variable time Ti the wavepackets partially overlap 
again and we apply a third pulse, of 6 fxs duration (a 7r/2 
pulse). This last pulse splits each wavepacket into the two 
momentum states. The interference of the overlapping 
wavepackets in each of the two momentum states allows 
the determination of the local phase difference between 
them. By changing the time T2 we vary Sx — xa—xb, the 
separation of ipA and tpB at the time of the final Bragg 
pulse. The set of data at different 5x constitutes a new 
type of spatial autocorrelation measurement that is sim- 
ilar to the "FROG" technique |17| used to measure the 
complete field of ultrafast laser pulses. From these mea- 
surements we obtain the phase profile of the wavepackets 
in the x direction. 




FIG. 2. (a-e) One of the two output ports of the interfer- 
ometer with To — 4 ms and Sx as indicated, (f) A plot of the 
density along the x direction of (d). 

Figure 0a-e shows one interferometer output port for 
different Sx (different T2) after an expansion time To = 4 
ms. In general, we observe straight, evenly spaced fringes 
(although for small Tq and T2 the fringes may be some- 
what curved). There is a value of Sx = xq ^ where 
we observe no fringes (Fig. ||c) and the fringe spacing de- 
creases as |<5x — xo| increases. Figure |^f, a cut through 
Fig. ||d, shows the high-contrast fringes |jl8). Our data 
analysis uses the average fringe period d, obtained from 
plots like Fig. ^f. 

The fringes come from two different effects: the inter- 
ference of two wavepackets with quadratic phase profile, 
and a relative velocity between the wavepackets' centers. 
The data can be understood by calculating the fringe 
spacing along x at output port 1 p9[. We assume that 
the phase cj> of the wavefunction /e^ can be written as 
<p = jx 2 + (3x. The equal spacing of the fringes im- 
plies, as predicted in the Thomas- Fermi limit |po| , that cf> 
has no significant higher-order terms |2JJ . The curvature 



coefficient a describes the mean-field expansion of the 
wavepackets and f3 describes a relative repulsion velocity. 
The velocity arises because the wavepackets experience 
a repulsive push as they first separate and again as they 
recombine. The density at port 1 (see Fig. 1) just after 
the final interferometer pulse is the interference pattern 
IVmi + ipBi\ 2 of the wavepackets ipAi and ipBV 

\f(x - Sx)e l ^ (x ~ Sx)2 - p{x - Sx)) + f(x)e L{ ^ x2+(3x) \ 2 , (1) 

where we assume that the amplitudes and curvatures of 
the wavepackets are equal and their velocities have equal 
magnitude and opposite direction. The cross term of (|l|) 
is 



2f(x — Sx)f(x)cos 



a Sx 



MSv 



x + C 



(2) 



where M is the sodium mass, MSv/h = 2/3, and C is 
independent of x [p2[ . Sv — vb — va is the relative re- 
pulsion velocity between the wavepackets ipAi and ipBi- 
Expression (|^) predicts fringes with spatial frequency, 



a Sx 



M Sv 



(3) 



where |k| = 2n/d. When there are no fringes, k = and 
the wavepacket separation Sx — xq = —MSv/ah. 
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FIG. 3. Plot of the spatial fringe frequency k versus Sx 
for To = 1 ms (filled circles) and 4 ms (open squares). The 
solid and dashed lines are linear fits to the data. 

Figure || plots the measured n vs. Sx Q for To = 
1 and 4 ms. The data are well fit by a straight line as 
expected from Eq. (3) in the approximation that a and 
Sv are independent of Sx. The slopes of the lines are the 
phase curvatures a, and the k intercepts give the relative 
velocities Sv. 

We checked the validity of the data analysis procedure 
by analyzing data simulated with a 1-D Gross-Pitaevskii 
(GP) treatment. Despite variations of Sv and a with 
Sx (due to their continued evolution during the variable 
time T2), we find that k is still linear in Sx. The slopes 
and intercepts in general are averages over the range of 
Sx used in the experiment. 

The interference fringes used to determine a and Sv are 
created at the time of the final interferometer pulse. Be- 
cause the two outputs overlap at that moment, we wait 
a time T3 for them to separate before imaging. Dur- 
ing this time, the wavepackets continue to expand. The 
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1-D simulations show that the fringe spacings and the 
wavepackets expand in the same proportion. We correct 
k (by typically 15 %) for this, using the calculated ex- 
pansion from a 3-D solution of the GP equation described 
below. 
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FIG. 4. (a) Plot of the phase curvature a versus the 
initial expansion time To showing the phase evolution from 
mean-field expansion. The solid line is a calculation using the 
Lagrangian Variational Method (LVM). (b) Plot of the rela- 
tive repulsion velocity Sv versus To. The solid curve is the 
calculated maximum repulsion velocity (when 5x = 0) and 
the dashed curve is the repulsion velocity averaged over the 
range of Sx used in the experiment. 

The different slopes and intercepts of the two lines in 
Fig. |3] show that the curvature a and relative velocity 
Sv of the wavepackets depend on the release time To be- 
fore the first interferometer pulse. Figure ||| plots the 
dependence of a and Sv on various release times To. The 
condensate initially has a uniform phase so that imme- 
diately after its release from the trap a = 0. We never- 
theless measure a nonzero a for To = ms because the 
BEC expands during T\ and T2. As a function of time, 
a behaves as D/D where D is the wavepacket diameter 
and D is its rate of change [Q . At early times when the 
mean-field energy is being converted to kinetic energy, D 
increases rapidly, increasing a. At late times, after the 
mean-field energy has been converted, D increases while 
D is nearly constant, decreasing a. 

We predict the time evolution of a using the La- 
grangian Variational Method (LVM) @. The LVM 
uses trial wavefunctions with time dependent parame- 
ters to provide approximate solutions of the 3-D time- 
dependent GP equation. In the model, the effect of the 
interferometer pulses is to replace the original wavepacket 
with a superposition of wavepackets having different mo- 
menta; e.g., the action of our first interferometer pulse is 
tpo - * (Vto + e t2kx ifj()) /V~2- We use Gaussian trial wave- 
functions in the LVM and, for simplicity, neglect the in- 
teraction between the wavepackets, to calculate the phase 
curvature a at the time of the last interferometer pulse. 
This result, with T\ = T2, is the solid line of Fig. 4a. 

We use energy conservation to calculate the rela- 
tive repulsion velocity Sv between ipAi and tpB 1 be- 
cause we neglect wavepacket interactions in the LVM. In 
the Thomas-Fermi approximation, we can calculate the 
amount of energy available for repulsion when Tq = 0. 



A trapped condensate has =// average total energy per 
particle, where n is the chemical potential After re- 
lease from the trap, it has |/z average mean- field energy 
per particle. Applying a ir/2 Bragg pulse to the BEC 
causes a density corrugation, which increases the mean- 
field energy to |/i per particle. In the approximation that 
the wavepackets do not deform as they separate and re- 
combine, one can show that 1/3 of the total mean-field 
energy goes into expansion of the wavepackets, and 2/3 
is available for kinetic energy of center-of-mass motion. 
Therefore jfj, of mean-field energy per particle is avail- 
able for repulsion. The corresponding repulsion velocity 
is only about 10~ 2 of a photon recoil velocity. The repul- 
sion energy and Sv decrease for larger To because both are 
inversely proportional to the condensate volume, which 
we calculate with the LVM. The two curves shown in Fig. 
||b are the calculated Sv when Sx = (solid curve) and 
Sv averaged over the different Sx used in the experiment 
(dashed curve). The 1-D GP simulations suggest that for 
small T , the results of the experiment should be closer 
to the solid curve; and for large T , closer to the dashed 
curve. The data is consistent with this trend. 
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FIG. 5. (a) Schematic representation of the interferometer 
in the trap, with the principle difference from Fig. 1 being the 
curved arrows indicating the acceleration of the wavepackets. 
(b) The relative velocity v between the two trapped wavepack- 
ets versus the interferometer time T. The solid line is a fit. 



In a related set of experiments we performed interfer- 
ometry in the trap. This differs from the experiments 
on a released BEC because there is no expansion before 
the first interferometer pulse j2f| and the magnetic trap 
changes the relative velocity of the wavepackets between 
the interferometer pulses (Fig. 5a). To better reveal the 
velocity differences, we choose T\ = T 2 = T to suppress 
fringes arising from the phase curvature. As with the 
released BEC measurements, we observe equally spaced 
fringes at the output of the interferometer, although the 
fringes are almost entirely due to a relative velocity v be- 
tween the wavepackets ipAi and ipm at the time of the 
third interferometer pulse. We obtain v from the fringe 
periodicity after a small correction for residual phase cur- 
vature po| . 

Two effects contribute to v. the mutual repulsion be- 
tween the wavepackets ipA and tps and the different ac- 
tion of the trapping potential on the two wavepackets in 
the interferometer. The latter effect occurs because after 
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the first Bragg pulse, ipA remains at the minimum of the 
magnetic potential while ipB is displaced. Wavepacket 
ipg therefore spends more time away from the center of 
the trap and experiences more acceleration than tpA- 

Following the last Bragg pulse, ipAl and 4>bi have a ve- 
locity difference which for our parameters can be approx- 
imated by v ps -^sin 2 (lu x T) + Sv Figure |b plots 
v versus T, and the curve is a fit to the above expression. 
We obtain the trap frequency uj x /2tt — 26.7(15) Hz, in 
excellent agreement with an independent measurement. 
We also obtain the relative velocity from the mean-field 
repulsion Sv = 0.49(12) mm/s, which we expect to be 
somewhat larger than for the released measurements be- 
cause the wavepackets contract, producing a larger mean 
field. 

In conclusion, we demonstrate an autocorrelating 
matter-wave interferometer and use it to study the evo- 
lution of a BEC phase profile by analyzing spatial images 
of interference patterns. We study how the phase curva- 
ture of the condensate develops in time and measure the 
repulsion velocity between two BEC wavepackets. Our 
interferometric method should be useful for characteriz- 
ing other interesting condensate phase profiles. For ex- 
ample, it can be applied to detect excitations of a BEC 
with characteristic phase patterns, such as vortices and 
solitons tljj2^ -|3l|] . The method should be useful for fur- 
ther studies of the interaction of coherent wavepackets 
and to study the coherence of atom lasers. 
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